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Abstract 

For some integrable systems, such as the open Toda molecule, the 
spectral curve of the Lax representation becomes the graph C = {(A, z) | 
z — A(X)} of a function A(X). Those integrable systems provide an in- 
teresting "toy model" of separation of variables. Examples of this type 
of integrable systems are presented along with generalizations for which 
A(X) lives on a cylinder, a torus or a Riemann surface of higher genus. 
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1 Introduction 

Separation of variables (SOV) is an old and new subject, which was born in the 
Hamilton- Jacobi theory of the nineteenth century p^, revived by Moser in the 
light of the modern theory of integrable systems [2], and extended by Sklyanin 
and his collaborators to a wide range of classical and quantum integrable sys- 
tems 0. Moreover, recent studies have revealed a number of new aspects of 
SOV in the context of bi-Hamiltonian structures 0J|S1|1>], universal Hamiltonian 
structures [7| , algebraic surfaces , geometric Langlands duality ^U] > string 
theory JT], affine Jacobi varieties |12M13 |. etc. 

This article is a supplement to the joint work with Takebe |2] on SOV of in- 
tegrable systems. This work was motivated by a rather naive question — what 
will be the simplest nontrivial examples of SOV of integrable systems? The 
first plan of this research was to use such systems as a "toy model" for better 
understanding of Sklyanin's work on SOV of quantum integrable systems (in 
particular, Calogero-Moser and Ruijsenaars-Schneider systems US]). Unfortu- 
nately, the program on quantum integrable systems has been unsuccessful due 
to technical difficulties. However, as we shall see below, searching for simplest 
examples of SOV is an interesting issue in itself. 
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2 Separation of variables of integrable Hamilto- 
nian systems 

Let us consider a Hamiltonian system with canonical variables q± , . . . , q n , px , . . . , p., 
and a Hamiltonian H = H(qx, . . . , q n ,Px, ■ ■ ■ ,p n )- This system is said to be 
integrable (or completely integrable) if it has n functionally independent (i.e., 
dHiA- ■ -AdH n ^ 0) first integrals Hi, ... , H n in involution (i.e., {Hj,Hk} = 0). 
The joint level set of the first integrals, 

Ht(q u ...,q n ,p 1 ,...,p n )=E i , £=l,...,n, (1) 

is a Lagrangian submanifold of the 2n-dimensional phase space, which (after 
canonical transformations of coordinates, if necessary) can be expressed as 

Pj = dS/dq, (2) 

by a function S = S(qi, . . . , q n , Ex, . . . , E n ). In the traditional terminology of 
the Hamilton-Jacobi theory, S is a "complete solution" of the joint Hamilton- 
Jacobi equations 

H i (qx,...,q n ,dS/dqx,...,dS/dq n )=E e , l=\,...,n. (3) 

The system can be thereby transferred to an action-angle systems with the 
action variables Eg and the angle variables 4>e — dS/dEp. 

In this setup, SOV means that the equations of the joint level set takes the 
"separated" form 

fj(\j,Hj,Ex,...,E n ) =0, j = l,...,n, (4) 

in another set of Darboux coordinates Xx, . . . , X n , (J-x, ■ ■ ■ , fJ-n ("separation vari- 
ables"). One can then assume the generating function S in the separated form 

n 

S = ^S,(A,, £!,...,£„) (5) 

i=i 

and obtain each term Sj by solving the ODE 

f j (X j ,dS j /dX j ,E 1 ,...,E n ) = 0. (6) 

Note that the existence of first integrals itself tells nothing about separability; 
a prescription of SOV has to be sought elsewhere. Sklyanin's "magic recipe" 
|3] is a prescription based on the Lax representation of integrable systems. The 
idea stems from the work of the Montreal group on generalized Moser systems 
[TC] , and is also closely related to the algebro-geometric method of the Moscow 
group |17| . Suppose that that the integrable system has a Lax representation 
of the form 

L(z) = [L(z),M(z)} (7) 
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with N x N matrices L(z) and M(z) that depends on a spectral parameter z. 
The characteristic polynomial 



N 

/(A, z) = det(AJ - L[z)) = X N +Y j n{z)X N - t (8) 

of L(z) does not depend on the time variable. A set of first integrals in involution 
can be obtained from the coefficients of expansion of pe(z)'s by a suitable basis 
of functions of z. The curve 

C = {(\,z)\f(\,z)=0} (9) 

defined by f(X,z) is called the "spectral curve" of the L-matrix L(z). At each 
point (A, z) of C is sitting the eigenspace of L(z) with eigenvalue A. For many 
integrable systems, this eigenspace is one-dimensional at almost all points of C, 
and the components of a suitably normalized eigenvector <^(A, z) are meromor- 
phic functions on C with a finite number of poles. It is the coordinates A,-, Zj of 
these poles that play the role of separation variables in Sklyanin's magic recipe. 
These variables are constrained by the equation 

f(Xj,Zj)=0 (10) 

of the spectral curve, and these equations can be interpreted as the separated 
equations (@J of the joint level set of first integrals. Moreover, in most cases, 
the Poisson brackets of these separation variables take the separated form 

{Xj, X k } = {zj,z k } = 0, {Xj, z k } = S jk g(X k ,z k ), (11) 

where <?(Afc, z k ) is a common function of Afe and z k (typically, 1, Afe, z k or X k z k ). 
In other words, the symplectic form of the system can be written as 

<12) 

3 If the spectral curve is the graph of a function 

Our concern lies in the case where the spectral curve becomes the graph of a 
function A(X): 

C={(X,z) \z = A(X)}. (13) 
Here is a list of examples of integrable systems of this type: 

1. Open Toda molecule 

2. Open Ruijsenaars-Toda molecule 

3. Rational and trigonometric Calogero-Moser systems 
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4. Rational and trigonometric Ruijsenaars-Schneider systems 

These integrable systems are usually understood to be associated with no spec- 
tral curve, because the usual Lax representation does not contain a spectral 
parameter. The fact is that they do have a Lax representation with a spectral 
parameter, and that the spectral curve takes the form mentioned above. 

SOV of these examples is almost parallel to more standard cases. The func- 
tion A(X) is expanded to a linear combination of a suitable basis /o(A), /i(A), . . . , /at (A) 
of functions as 

N 



A(A) = /o(A) + 5^/,(A), (14) 
1=1 

and the separation variables comprise N pairs (A 1; z±), . . . , (Ajv, of variables 
that satisfy the equations 

z j = A(X j ), j=l,...,N. (15) 

These equations, in turn, determines the first integrals Hg as a function Ei 
II, ; A ■ . . . . , Aat, zi, . . . , zn) of the separation variables. More precisely, we as- 
sume that the non-degeneracy condition 

A = det(/ fc (A i ) \ j,k=l,...,N)^0 (16) 

is satisfied. One can then use Cramer's formula to solve the foregoing equations 
for the coefficients Ei as 

Et = A t /A, (17) 

where Ai denotes the determinant in which the £-th row of A is replaced as 
ae(Xj) — > Zj — ao(Xj). This gives an explicit formula of the first integrals. 

For illustration, we consider the case of the open Toda molecule in detail 
below. 



4 Example: open Toda molecule 

4.1 L-matrix and spectral curve 



The open Toda molecule of length N has the Hamiltonian 

N N 

2 



N N 

^^E^+E^- 1 ^- ( 18 ) 

j'=i i=2 

A Lax representation with a spectral parameter z is presented by Krichever and 
Vaninsky ^H] . The L- matrix reads 

N N N 

L{z) = X c o E o,3-i + E M-',, + E E i-U + zEn ^ 

3=2 , I 3=2 
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where Ej^ stands for the matrix unit with the (J, k) element equal to 1 and 
the others equal to 0; bj and Cj are the Flaschka variables, bj = pj and Cj = 
c ij-i-<ii , Note that L = L(0) is the usual tri-diagonal L-matrix without spectral 
parameter. 

This _L-matrix can be derived from the L- matrix L c \(z) of the closed Toda 
molecule (namely, the Toda lattice with periodic boundary condition qj+N = Qj, 
Pj+N = Pj) as follows. The Hamiltonian of the closed Toda molecule is the sum 

H cl = H + c u ci-e""-" 1 (20) 

of the Hamiltonian of the open Toda molecule and the potential between qjy 
and qi. The associated L-matrix is also the sum of two pieces: 

L c i(z) =L(z) + c 1 z- 1 E 1 m. (21) 

In the limit as the coupling between q^ and q\ is turned off (c\ — > 0), this 
L-matrix becomes the aforementioned matrix L{z). 

This correspondence carries over to the spectral curve as well. The charac- 
teristic polynomial of L c \(z) takes the well known form 

det(A7 - L cl (zj) = z + a ■ ■ ■ c N z- x - P(X), (22) 

where P(A) is a polynomial of degree N. The spectral curve defined by this 
polynomial is a hyperelliptic curve of genus N — 1. In the limit as c\ — > 0, the 
characteristic polynomial reduces to 

det(A7 — L(z)) = z — A(X), (23) 

where A(X) is nothing but the characteristic polynomial of L = L(0): 

A(X) = det(A7 - L). (24) 

The spectral curve of the open Toda molecule thus turns out to be the graph of 
A{X). 

4.2 Separation variables 

Having the i-matrix L(z), one can obtain separation variables from poles of a 
suitably normalized eigenvector of L(z). This is indeed achieved in the work of 
Krichever and Vaninsky 18, . One can derive the same set of separation variables 
from those of the periodic Toda molecule; all necessary data for the latter can 
be found in the work of Flaschka and McLaughlin ^5]. We here present yet 
another approach based on Moser's classical method for solving the open Toda 
molecule [2*D], 

Moser's method uses one of the diagonal elements of the resolvent (XI— L) : 
R(X) = t e N (XI - L)~ 1 ejy, *e N = (0, . . . , 0, 1). (25) 
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This is a rational function of the form 

R(X) = B(X)/A(X). (26) 

Since the numerator B(X) is a polynomial of degree N—l, R(X) can be expanded 
to partial fractions as 

N N 

i? ( A )=El^' A(X) = l[(X-a J ). (27) 

3 = 1 3 j=l 

Moser discovered that the new dynamical variables aj and pj obey the simple 
evolution equations 

aj=0, Pj = a jPj , (28) 

which can be readily solved as a,- (*) = aj(0), pj(t) = e^^pj. The L-matrix 
can be reproduced from ctj and pj by the Stieltjes method of continued fraction. 
This is the way how Moser solved the open Toda molecule. 

From our point of view, the roots of B(X) are a half of separation variables: 

N-l 

B(X)=l[(X-X k ). (29) 

k=l 

The other half z\, . . . , zjv-i are given by the values z k = A(Xk) of A(X) at the 
roots of B(X). The Poisson brackets of these variables turn out to take the form 

{Xj,Xk} = {zj,Zk} = 0, {Xj,Zk} = SjkZk- (30) 

One can derive these Poisson brackets in several different ways. The most 
direct way is to use the result of Faybusovich and Gekhtman j^] on the Poisson 
brackets of A(X) and B(X). 

A comment on the number of variables will be in order. These 2N — 2 
separation variables Xf., Zk amount to N — 1 of the N degrees of freedom of the 
open Toda molecule. The rest is the degree of freedom carried by the center of 
mass Q — J2f=i Qj/N and the total momentum P = J^f^iPji which disappear 
in the center-of-mass frame. 

4.3 Many faces of symplectic form 

We now have at least three different sets of dynamical variables besides the 
canonical ones qj,Pj — Moser's variables <x,-, pj, the separation variables Xk, Zk, 
and the action-angle variables Ei,(j)£. Accordingly, the symplectic form £1 of the 
system has many different expressions. We present their explicit form below (in 
the center-of-mass frame, for simplicity). 
Let us start from the expression 

N-l 

n = dlog A(A fe ) A dX k (31) 
fc=i 
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in the separation variables and Zk = B(\k). Substituting A(Xk) — Yif^ii^k — 
aj) leads to an intermediate expression 

k=l j=l ■> 

which can be reorganized to the expression 

N 

CI = ^ d\ogB(aj) A daj (33) 

in the new variables aj and B(aj). As Faybusovich and Gekhtman point out |21| 
(see also Vaninsky's recent paper[22]), this is exactly the symplectic form that 
Atiyah and Hitchin introduced on the moduli space of SU(2) monopoles. In 
terms of Moser's variables aj and Pj = B(aj)/A'(ctj), the last expression reads 

= f>log Pj A da, + J2 dajAdak . (34) 

i - / 1 ^7 ^k 

j=l j=ik J 

An expression in action-angle variables can be obtained as follows. Let En 
be the coefficients of A(X); 

N 

A(\) = \ N + Y,EiX N ~ e . (35) 

1=2 

[Ei is equal to the total momentum, hence absent in the center-of-mass frame.) 
The foregoing expression of f2 can now be rewritten as 

N 

n = ^2dE e Adfa, (36) 

i=2 

where 



£/ m iX - (37) 

This shows that fa is an angle variable conjugate to the action variable Eg. It 
is easy to see that generating function S of the action-angle variables, too, has 
a similar simple expression: 

S=J2 \og A(X)d\. (38) 
fc=i ^ 
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5 Other examples 



5.1 Open RT molecule 

The open Ruijsenaars-Toda (RT) molecule may be thought of as a limit of the 
closed RT molecule (namely, the RT lattice with periodic boundary condition 
1j+N = Qj and Pj+n = Pj)- The passage to the open molecule is parallel to the 
case of the Toda lattice. The closed molecule has the Hamiltonian 

N 

H cl = J2 eP3 ( l + e ' l3 ~ 1 ~ q3 ) 1/2 ( l + e«^+ 1 ) 1/2 , (39) 
i=i 

which turns into the Hamiltonian H of the open molecule as the potential e qN ~ qi 
between q\ and qw is turned off. 

SOV of this case is more delicate than the Toda molecules, because curves 
arising here are NOT a spectral curve in the usual sense. This is already the 
case for the closed molecule. The closed RT molecule was solved by Bruschi and 
Ragnisco |24| by an algebro-geometric method, but they had to to use an unusual 
Lax representation with a matrix L(X, z) that depends on two parameters A and 
Z. Since L(X, z) amounts to L(z) — XI in the usual Lax representation, they 
considered the curve defined by the equation 

f(X,z) = det L(X,z) = 0, (40) 

and obtained a set of separation variables as Flaschka and McLaughlin ^] did 
for the closed Toda molecule. Therefore, bearing in mind the case of the Toda 
molecules, it is not difficult to see what appears in the limit to the open RT 
molecule. 

In the limit to the open molecule, as expected, the hyperelliptic curve of 
Bruschi and Ragnisco degenerates to the graph of a polynomial A(X) of the 
form 

JV-l 

A(X) = X N + E t X N - 1 + 1. (41) 
i=\ 

Separation variables Ai, . . . , Ajv-i are the roots of a polynomial B(X) of degree 
N — 1, and their conjugate variables z%, . . . , zjv-i are defined by Zk = A(Xk)- 
The Poisson brackets are slightly different from those of the Toda molecules: 

{Xj, A fe } = {zj, z k } = 0, {Xj, z k } = S jk X k z k . (42) 

The symplectic form Q of the system, too, is modified as 

N-l 

Q = d\ogA(X k ) Ad\ogX k 

fc=i 

N 

= ^ d log B(Xj ) A d log «j 
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N-l 



= j2 dEeA ( 43 ) 



1=1 

where ay's are the roots of A(A), and <fits are given by 



5.2 Rational and trigonometric CM systems 

The rational and trigonometric Calogero-Moser (CM) systems may be thought 
of as degeneration of the elliptic CM system. This allows us to find a Lax 
representation with a spectral parameter for these systems. 

Let us recall the following L-matrix for the elliptic CM system: 

N 

L(z) = ^2pjE jd + giH<lj - Ik, z)E j}k , (45) 

j=l j^ik 

where 4>{u, z) is written in terms of the Weierstrass a function as 

a(u + z) , 
4>(u,z)= . 46 

<7(U)(J{Z) 

As the elliptic curve degenerates to a rational curve, (j)(u, z) turns into a trigono- 
metric (hyperbolic) function 

sinh(u + z) 

(f>{u, z) — » — r~r~\ n~r — coth(w) + coth(z) 

smh(tt) smh(z) 

or a rational function 

u + z 1 1 



<t>(u,z) 

Remakably, the limit of <j)(u, z) splits into two functions of a single variable. 
Accordingly, the limit of L(z) becomes the sum of two matrices as 

, > J L + gicoih(z)K (trigonometric) , 
L(Z > ^ { L + giz- Y K (rational) (47j 

where 

N 

L = J2 Pj E 3^ + X! gi cotn (<2j - 1k)Ej,k (48) 

3 = 1 j^k 

in the trigonometric limit, 

N 

L = J2 p i E ™ + X! • //: ' / < - 9h)~ lE i,h (49) 

3=1 j^k 



9 



in the rational limit, and 

K = Y,E j , k (50) 

for both cases. L is nothing but the usual L-matrix of the trigonometric and 
rational CM systems; K is the constant matrix that appears in the method of 
Hamiltonian reduction (or the "projection method") for solving these systems 
|2(il I27| . The linear combination of L and K in the limit (|47|l of L(z) gives 
an L-matrix with a spectral parameter for the rational and trigonometric CM 
systems 

We modify these rational and trigonometric L-matrices slightly. Note that 
K can be expressed as 

K = e t e-I, *e= (1,...,1). (51) 

Therefore the term proportional to K in l|47ll is a linear combination of a rank- 
one matrix and a scalar matrix; the latter does not affect to the Lax equation. 
We can thus drop the scalar term and redefine the L-matrix for the degenerate 
CM systems as 

L(z) = L + gicoth(z)e t e (52) 
for the trigonometric case and 

L(z) = L + giz~ 1 e t e (53) 

for the rational case. 

We now apply the so called Weinstein-Aronszajn formula 

det(M + x*y) = (1 + V-M" -1 *) det M 
to the characteristic polynomial of L(z). This yields the expression 

det(A7 - L(z)) = { " gl cot h(^i( A ) (trigonometric) 

I "l(Aj — giz "i(A) (rational) 

where 

P (A) = det(A7 - L), P X (A) = t e(A/^i)e. (55) 

(M denote the cofactor matrix of M.) The equation of the spectral curve of 
L(z) thus turns out to take the form 

for the rational case and 



,,22 



e 



P (A)+Pi(A) 
Po(A)-Pi(A) 



(57) 
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for the trigonometric form. This is what we have sought. 

It should be mentioned that the same equation of the spectral curve for the 
trigonometric case was also derived by Vaninsky from a different point of 
view. We can proceed further to the construction of separation coordinates, 
though we omit details here. 



5.3 Rational and trigonometric RS systems 

The rational and trigonometric Ruijsenaars-Schneider (RS) systems can be 
treated in the same way as the CM systems. We again start from the elliptic 
RS system, and consider the limit as the elliptic curve degenerates to a rational 
curve. It is convenient to use the following L-matrix of the elliptic RS system 
(which is gauge equivalent to the L-matrix of Bruschi and Calogero |29j): 

N 

L(z) = {hjhk) 1 ' 2 ^ - q k + 7, z)Ej, k , (58) 

3,k=l 

where 7 is a constant and hj 's are fairly complicated functions of the canonical 
variables. In the limit to the trigonometric or rational models, the L-matrix 
again splits into the sum of two matrices, one being the L-matrix L without 
spectral parameter, and the other a rank-one matrix. One can thus use the 
Weinstein-Aronszajn formula again to show that the spectral curve is the graph 
of a function A(X). The same curve was derived by Braden and Marshakov [5U| 
in the context of supersymmetric gauge theories. 



6 Variants and generalizations 

The goal of this research is to search for new integrable systems from the point 
of view of SOV. Let us show a few results below. 



6.1 Trigonometric and elliptic analogues 

Apart from the CM and RS systems, the foregoing examples are associated with 
a pair of polynomials A(X) and B(X). The main result of the joint work with 
Takebe [2] is to construct the following analogues defined on a cylinder and a 
torus: 

1. Trigonometric version 

N N 

A(X) = J[ sinh(A - aj ), B(X) = JJ sinh(A - Xj), (59) 

2. Elliptic version 

N N 

A(X) = l[ ( r(X-a j ), B(X) = JJ a(X - Xj). (60) 
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The construction can be further extended to the case where A(X) and B(X) have 
a different number of zeroes. 



6.2 Elliptic fibration 

Another possible elliptic analogue is to modify the z direction of the (A, z) space 
to an elliptic curve. A naive way will be to replace 

N N dz 

= d log Zj A dXj -> = — - A dXj (61) 

3=1 i=i "' 

where yj and Zj are constrained by the equation 

Vj = 4z| + 52 + .93 

of an elliptic curve. Actually, one can go further by allowing the coefficients gi 
and <73 to depend on Xj as 

y] = 4zf + g 2 (X j )z j +g 3 (X j ), (62) 

where 32(A) and 33(A) are functions (e.g., polynomials) of A. Geometrically, this 
is the structure called an "elliptic fibration". Integrable systems of this type, 
too, are considered in the joint work with Takebe |14| . 



6.3 Higher genera 

Generalization to a complex algebraic curve Co of genus g > 1 (and a Riemann 
surface of more general type) will be achieved in several different ways. Let us 
propose an idea below. 

The idea is to choose the (A, B) pair as nonzero elements of the vector 
space C{R\ + • • • + Rn) of mcromorphic functions f(P) (P € Co) on Co with 
at most simple poles at N fixed points Rx, . . . , Rjy. If N > 2g — 2 (or if 
N > g and Ri, . . . , Rpj are in general position), the Riemann- Roch theorem 
shows that this vector space is spanned by N — g + 1 linearly independent 
elements / (P), A(P), . . . , f N - g (P)- We write A(P) and B(P) by this basis as 

N-g N-g 
A(P) = Ulh{P), B(P) - ]T Vifi(P) (63) 

1=0 i=0 

and view {ui)^~q and (ve)fS 9 as homogeneous coordinates on two copies of 
the projective space P"" 9 . The (A, B) pair is thus parametrized by 2A^ — 2g 
"moduli" , which play the role of dynamical variables. Now let us choose a 
holomorphic differential dX on Co, and consider its (multi- valued) primitive 
function 

A(P) = / dX. (64) 
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In this quite general setup, one can use a trick based on the residue theorem 
[T4*] to derive the identity 



N N 



J2dlogB(P 3 ) AdX(Pj) = ^d\ogA(Qj) AdX(Qj), 



(65) 




where Pj and Qj are zeroes of A(P) and B(P), respectively; d is understood to 
be the total differential in the 2N — 2g moduli of the (A, B) pair. This defines 
a symplectic form f2 on the moduli space. The inhomogeneous coordinates 
Hi — Ui/uo give a set of functions in involution with respect to f2, thereby 
defining an integrable system. The role of separation variables is played by 



though they are not independent variables. 
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